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Abstract 

We suggest that a singlet fermion S exists beyond the standard see-saw struc- 
ture. It mixes with light neutrinos via interactions with the right-handed neutrino 
components, so that v e — > S conversion solves the solar neutrino problem. Su- 
persymmetry endowed with R-symmetry is shown to give a natural framework for 
existence, mass scale and mixing (sin 2 2# es ~ (0.1 — 1.5) ■ 10~ 2 ) of such a fermion. 
Models with an approximate horizontal symmetry are constructed, which embed 
the fermion S and explain simultaneously solar, atmospheric, hot dark matter 
problems as well as may predict the oscillation —>■ v e in the region of sensitivity 
of KARMEN and LSND experiments. 



1 Introduction 



The solar neutrino problem Q, the deficit of muon neutrinos in atmospheric neutrino flux 
||, the large scale structure of the Universe || and possible candidate events in a search for 
z/ M — > v e oscillations [|J (see however H) give indications on non-zero neutrino masses and 
lepton mixing. Simultaneous explanation of all (or some) of these problems may call for the 
existence of more than three light neutrinos which mix among themselves ||. Strong bounds 
on the number of neutrino species both from the invisible Z°-width and from primordial 
nucleosynthesis (NS) (7| require the additional neutrino to be sterile (singlets of SU (2) x U{1)). 
The right-handed (RH) components of known neutrinos are natural candidates for such sterile 
states. However, in such a case one has to depart from the conventional see-saw mechanism 
which implies large masses to the RH components. 



A number of schemes with light sterile neutrinos has been suggested |p|- ||12|| . Most of 
them are based on radiative mechanism of mass generation or on some hybrid schemes which 
include both the elements of the see-saw and radiative mechanisms. In these schemes sterile 
neutrino is considered on the same footing as the usual neutrinos. The lepton number is 
broken typically at the electroweak scale. 

We will consider another possibility. We suggest that usual see-saw mechanism works with 
all three right-handed neutrinos having large Majorana masses: M ~ 10 10 — 10 12 GeV. At the 
same time the theory contains an additional singlet fermion S which has its origin beyond the 
standard lepton structure. The singlet S is very light and mixes with neutrinos. 

Supersymmetry (SUSY) can provide a natural justification for the existence of S. Many 
extensions of the standard model contain singlet scalar fields: singlet majoron ||13|| , invisible 
axion JT4|, or scalars for spontaneous generation of the /i-term [15], etc.. The supersymmetric 
partners of such scalars could be identified with S. Moreover, SUSY can play a crucial role in 
the determination of mass scales in the singlet sector. 

In this paper we consider possible origin of light fermion S, its mass and mixing with light 
neutrinos. The models with S are constructed so that they can simultaneously explain the 
above mentioned neutrino anomalies. 
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2 Light singlet fermion and the solar neutrino problem 



Primordial nucleosynthesis (as well as the data from SN87A) gives strong bound on the oscilla- 



tion of active neutrinos into sterile neutrino [IB|. This practically excludes v^ — > S oscillations 
as a solution of the atmospheric neutrino problem. The singlet fermion with mass in the 
eV-range could be considered as a candidate for hot dark matter (HDM) f6j. However, if its 
density satisfies the NS bound on the number of additional neutrino species: 5N U ~ 0.1, it 
can not reproduce the optimal parameters || for the large scale structure formation in the 
Universe: m$ ~ (2 — 5) eV and Q s ~ 0.2, where fl s is the energy density of S in the Universe 
in the unit of the critical density. Therefore it may happen that the only place where singlet 
fermion plays a role is the solar neutrino problem. 

Let us find the region of parameters Am 2 and sin 2 29 es , where 9 es is the mixing angle of 
v e with S for which the resonance conversion v e — > S inside the Sun can explain the existing 
data. It is instructive to compare the sterile {y e — > S) and the active neutrino (v e — > Uf) cases. 
The v e — > S solution of the solar neutrino problem differs from the v e — ■> v^ (v T ) solution in 
two ways. 

(1) The effective density, p s , for v e — > S conversion is smaller than that, pf, for v e 



conversion: 

n Y — -Y 

l J s 1 e 2 n 

7 f = Y e ■ 

Here Y e and Y n are the number densities of the electron and the neutron per nucleon, respec- 
tively. In the center of the Sun one gets p c s / p°f — 0.76 JT?]. The central density p c determines 
position of the adiabatic edge of the suppression pit: (E/Am 2 ) a oc l/p c . Consequently, in the 
v e — > S case the adiabatic edge is shifted to larger E/Am 2 in comparison with the flavour 
case: (E/Am 2 ) s = (E / Am 2 ) ^ ■ pj/p s . The position of the nonadiabatic edge depends on p/p 
and the difference between the flavour and sterile cases is practically negligible. 

In the region of small mixing solutions [], the allowed values of Am 2 are determined essen- 
tially by (E/Am 2 ) a and by Gallium experiment data |18|]. Therefore the shift of the adiabatic 
edge for v e — > S to larger E/Am 2 results in corresponding shift of Am 2 to smaller values: 



Am 2 | s ~ ^ Am 2 \ f ~ 0.76 Am 2 \ f . (1) 
P C f 

(2) The fermion S has no weak interactions, and therefore S flux from u e —* S conversion 
does not contribute to the Kamiokande signal {ve ve scattering) in contrast with flavour 
1 Large mixing domain is excluded by primordial nucleosynthesis data 
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case, where interacts via neutral currents. This influences the allowed region of mixing 

angles. Indeed, for unfixed original Boron neutrino flux (which has the largest, ~ 50%, 

theoretical uncertainties) the bound on sin 2 29 es is determined by the "double ratio" [TS | : 

_ Rat 
tiHiK = -5— , 

where Rat = Q°Ar/Qjw M an d R» e = $°b s /$b SM are the suppressions of signals in Cl-Ar 
and Kamiokande experiments, respectively. Here Q^f M j ^ B SM are the predictions in the 
reference model (e.g. (Tj|) and Q° A b ^ , <&£ S are the observable signals. Due to the (z/ T )-effect, 
Rf, in the sterile case is smaller than R[ in the flavour case, and since Ra t is the same in 
both cases one gets R S H / K > Rh/k- With diminishing of 9 es , the suppression of due to 
conversion weakens and the effect of (z/ T ) decreases, therefore R S H / K approaches R h /k- As 
a consequence, the lower bound on 9 es coincides practically with that for flavour conversion: 
sin 2 2# e <i ~ (0.8 — 1.0)-10~ 3 [18]. On the contrary, with increase of 9 es the suppression of $b due 



to the conversion becomes stronger, so that Rat — ► and R s Ue — > 0; at the same time due to 
the neutral current effect of v^{v T ), Rf e - 0.16. Therefore R f H/K - 0, whereas R' H/K does not 
change strongly. We have found R S H/K ~ 0.77, 0.74, 0.72, 0.69 for sin 2 29 es ~ 2 • 10" 3 , 5 • 10~ 3 , 
10~ 2 , 2 • 10~ 2 , respectively. The experimental value of the double ratio is Rh/k = 0.67 ±0.11. 
However for large sin 2 29 es the original flux of Boron neutrinos should be large (to compensate 
for strong suppression effect). If we restrict < 1.5$fj SM , then the bound on the mixing 
angle becomes: sin 2 2# es ~ 1.5 • 10 -2 . This also satisfies the NS bound Wfi . 



Resonance conversion implies that mg > m Ue and if there is no fine-tunning of masses, 
nis — \J Am 2 . Thus using (P and known results for flavour conversion as well as bounds on 
sin 2 29 es discussed above we get the following range of the parameters: 

m s ~ (2-3)-10" 3 eV 
sin^ es ~ tan^ es ~ (2 - 6) ■ 10" 2 . (2) 

3 Mass and mixing of singlet fermion via right-handed neutrino 

Let us consider the following Lagrangian, 

m M 
C = j^L e u c e H 2 + -^vXe + ™ es v c e S , (3) 

where L e is the lepton doublet, H 2 is the Higgs doublet and v c e is the right-handed neutrino 
component. We suggest that there is no direct coupling of S with L e due to a certain symmetry, 
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and the mass term SS is absent or negligibly small. The Dirac mass m e and the mixing mass 
m es are much smaller than the Majorana mass M e : m e , m es << M e . The Lagrangian (||]) leads 
to the mass matrix in the basis (S, u e , v^)\ 

I m es \ 
M= m e (4) 
V m es m e M e J 

The diagonalization of (^j) is straightforward: one combination of the v e and S, 

uq = cos es v e + sin 9 es S , 
is massless, and the orthogonal combination, 

Vi = cos 6 es S — sin 9 es v e , 
acquires a mass via the see-saw mechanism: 

mi jg^ ■ (5) 

The mass of the heavy neutrino is ~ M e . The v e -S mixing angle is determined by 

tan6U = — - , (6) 

m es 

and correspondingly sin 2 29 es = 4:[m es me/ml s + ml} 2 . Taking for m e the typical Dirac mass of 
the first generation: m e ~ (1 — 5) MeV, and suggesting that v e — > 5* conversion explains the 
solar neutrino problem with mi = m$ as in (^|), we find 

m es = ~ (0.02 - 0.3) GeV . (7) 

tan v es 

According to (|) the RH mass scale is 



77? 

M e ~ m 2 /mi = - (10 8 - 3 • 10 10 ) GeV . 

mi tan 9 e 
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'es 



Consider now the models which lead to the Lagrangian ([|) with parameters (0) and (§). 
The simplest possibility is to use the U{1) symmetry of lepton number and to generate 
the masses in by VEV (a) of the scalar singlet, a. Prescription of the lepton charges 
(1, —1, —3, 2) for (z/ e , v% S, a) admits the following interactions in the singlet sector: 



o 2 ,„_a 3 



C = hvlvla + h'vtS— + h"SS^ , (9) 
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where Mp\ is the Planck mass. The Lagrangian @ reproduces the mass terms of (|3|) with 
M e = h(a) and m es ~ h! (a) 2 jMp\. The desired values of M e and m es (0) @ can be achieved 
with e.g., (a) ~ 10 10 GeV, h ~ 1 and /i' ~ 10~ 2 . The last term in (Q) generates the Majorana 
mass of S, mss = h" (a) 3 /Mpi, so that all the neutrinos are massive. For h" ~ one gets 
m\ ~ m 2 es /M e as before, whereas the smallest mass is ~ m^ssML < m i- 

Let us consider the possible role of supersymmetry in the appearance of the singlet fermion 
and in the determination of its properties. In principle, S can be a superpartner of the gold- 
stone boson which appears as a result of spontaneous violation of a certain global symmetry 
like lepton number or Peccei-Quinn symmetry. In this connection, let us consider a SUSY 
model with spontaneous violation of lepton number. The superpotential of the singlet ma- 
joron model is 

W = jrfrLtfHi + jv%vlo - X(aa' - M 2 )y , (10) 

where lepton numbers of the superfields (L e , z/;, a, a', H 2 , y) are (1,-1,2,-2,0,0). Lepton 
number is spontaneously broken by non-zero VEV's of a and a'. As the result, the majoron 
and its fermionic partner, the majorino, are massless in the supersymmetric limit. 

The identification of the majorino with 5* requires, however, the following complication of 
the model. 

(1) Supersymmetry breaking results in appearance of non-zero VEV of y which generates 
the mass of the majorino S: mss — Mv)- The soft-breaking terms X(A y aa' — B y M 2 )y + h.c, 
where A yi B y ~ 0(7713/2) are soft-breaking parameters give 

(y)^^(A y -B y ), (11) 

and consequently too big value of m S s — (A y — B y )/2 ~ 0(m 3 / 2 ), whereas A y — B y ~ 10~ 3 
eV is needed. One can get A y — B y = at tree level in no-scale supergravity or in the case 
of the non-minimal kinetic term discussed below. However, non-zero value of A y — B y will be 
generated due to renormalization group evolution of soft-terms. In order to suppress the mass 
below the solar neutrino mass scale (0), tuning of parameters is needed: / ~ 10~ 5 , if all three 
generation of leptons are taken into account. 

(2) Mixing of neutrinos with the majorino implies violation of R-parity. In (|T0"| ) the mixing 
can be induced by the second term if sneutrino v c gets non-zero VEV {vf). The latter requires 
the introduction of terms like v c e F[Xi) + W'(Xi), where F(Xi) and W'(Xi) are the functions 
of new superfields Xi. They should be arranged in such a way that in the global SUSY limit 
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F does not get a VEV: (F(Xi)) = 0, and after lepton number breaking linear term ~ M 2 v c e 
appears in the superpotential. Then the corresponding soft-term will give the mixing mass. 
One can find that the additional sector requires at least several new fields with non-zero lepton 
numbers which leads to further complication. R-parity violation is a general feature of models 
in which S is identified with fermionic superpartner of scalars acquiring non-zero VEV as in 
models for majoron, axion and /i-term. 

The above problems can be avoided in models with R-parity conservation. In this case, 
the lightest supersymmetric particle can be served as cold dark matter of the Universe. To 
preserve R-parity one should place the singlet S in the superfield with zero VEV. Consider 
the superpotential: 

777 A 

W = jT^L e u c e H 2 + jvlvla + fvlSy - -{a 2 - M 2 )y . (12) 

Its structure is determined by the R-symmetry under which the fields carry the R-charges: 

(1,1,-1,2,0,0) for (L e ,v c e ,S,y,a,H 2 ) . 

Note that the R-symmetry forbids the bare mass terms SS as well as the coupling SSa. Since 
lepton symmetry is explicitly broken no majoron appears. In the global SUSY limit, a gets 
non-zero VEV (a) ~ M ~ 10 11 GeV which generates the Majorana mass of v c e : M e = f{a). 
SUSY breaking induces the following soft-breaking terms in the scalar potential: 

nrr) \ 

Vsoft = {A i ^L e ^ 2 +/A^>>+fA 5 ^^--(^a 2 -5,M 2 )y+h.c.}+^m i 2 |^| 2 , (13) 

where z t denotes the fields appearing in the superpotential (|T2|) and A L , etc., are the soft- 
breaking parameters. Minimization of the potential shows the following: 

(1) The fields L e , v c el S do not develop VEV and therefore R-parity is unbroken. 

(2) The field y acquires non-zero VEV due to the soft-breaking terms as in (pTTj) . Consequently, 
the mixing mass for S and v c e appears: 

m es = ^(A y -B y ) (14) 

Since m es » m 1; no strong tunning of A y — B y is needed as in the previous case (PH). 
At A y — By ~ 0(777.3/2), the desired value of m es (0) can be obtained by choosing f'/X ~ 
10~ 3 — 10~ 2 . However, more elegant possibility is that A y = B y at the Planck scale but a 
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non-zero value for A y — B y is generated due to renormalization group evolution through the 
differences in interactions of a and y. In this case one expects 

A 2 , . 

where A represents a combination of the constants A, / and /'. As a consequence, the value 
m es ~ 0.1 GeV does not require smallness of A or /'. 

The equality A y = B y at the Planck scale can be achieved by the introduction of non- 
minimal kinetic term with mixings between the observable and hidden sectors. Let us introduce 
the following Kahler potential: 

K = CC + CC{a^- + a-f-) + ZZ , (16) 

Mpi Mpi 

where C and Z represent an observable and hidden sector field, respectively. Then usual 
assumption that the observable sector has no direct coupling to the hidden sector in superpo- 
tential, W = W(C) + W(Z), leads to the universal soft-terms: 

V soft ~m 3/2 W(C)+h.c. , (17) 

provided a = (W(Z)) / (M Pl dW/dZ + W(Z)Z/M Pl ). Note also that the field C does not 
acquire a soft-breaking mass. This mechanism can be generalized to arbitrary number of 
observable sector fileds. For our purpose C = a,y, i.e., we couple a and y to the hidden sector 
field Z with the above-mentioned choice for a. 

Note that a field plays two-fold role in the model: it gives Majorana mass of v c and it also 
generates mixing of v c with S by inducing a VEV for y after the SUSY breakdown. Moreover, 
a can be used to generate the //-term via the non-renormalizable interaction: 

2 

n H X H 2 . (18) 



M Pl 

The //-term can also be generated through the renormalizable interaction: yH\H2 in the case 
of (y) ~ C(m 3/2 ). 

It is easy to incorporate the spontaneous violation of lepton number or/and Peccei-Quinn 
symmetry into the model. As in (|T0| ) one should introduce the superfield a' with lepton number 
—2 and zero R-charge and replace the a 2 term of (0) by oa' . In this way the //-term flT8|) can 
be naturally related to the solution of the strong-CP problem via Peccei-Quinn mechanism 
21 1, and the majoron will coincide with the invisible axion 
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4 Models with light singlet fermion 

Two other neutrinos, and u T , can be included in the scheme by adding to (HD analogous 
terms with L a and v c a (a = [i, r). Then as in (|]), mixing of these neutrinos with singlet S is 
determined by tan(9 as = m as /m a , where m as and m a are the corresponding mixing and Dirac 
masses. 

Primordial nucleosynthesis gives strong bounds on the angles 9 as and/or on masses of light 
neutrino components: ~ m 2 a /M a . Suppose that S is family blind, and its couplings with all 
neutrinos are universal: m es ~ m^ s ~ m TS ~ (0.02 — 0.3) GeV. Note that this mass scale 
(motivated by solar neutrinos) is of the order of Dirac masses in the second generation. Then 
with m^ s ~ m M ~ 0.3 GeV, one gets tan 6^ s ~ 1, and if m 2 — rri^/M^ ~ 0.1 eV, the oscillation 

— > S could explain the deficit of atmospheric neutrinos. However, this possibility is strongly 
disfavoured by NS data. For m fl ~ 1 GeV one has sin 2 29^ ~ (0.2 - 4) ■ 10~ 2 , and the NS 
bound [0] is satisfied if Am 2 ~ (10~ 4 — 10~ 3 ) eV 2 , or m 2 ~ 3 • 10~ 2 eV. For the third generation 
(m T ~ 100 GeV), analogous figures are: sin 2 26 TS ~ (0.2 — 5) • 10~ 6 and m 3 < 3 eV. Therefore, 
the cosmologically interesting masses of v T are admitted. Note that the bound on 777,2 form NS 
and values of m\ and 7773 desired by solar and HDM problems can be reproduces by moderate 
mass hierarchy of the RH neutrinos: M a ~ 10 10 — 10 12 GeV. 

To have simultaneously neutrinos as HDM and the solution of the atmospheric neutrino 
problem via — > v T oscillations one needs — m 3 ~ 2 eV. In this case (Am 2 ~ 4eV 2 ) the 
NS bound: sin 2 29^ s ~ 10 -6 implies m^ s /m^ < 5 • 10 -4 or m^ s ~ 0.5 MeV at ~ 1 GeV, i.e. 
the coupling of S with u e should dominate: 777 es >> m^. 

Both the dominance of S—v° coupling and the near degeneracy of neutrinos corresponding 
to the second and the third generations {m<i ~ 7773) can arise as consequences of some family 
(horizontal) symmetry. 

Let us consider [7(1 ^-symmetry with charge prescription (0, 1, —1) for the first, the sec- 
ond and third generations of leptons, respectively. Each generation includes the left-handed 
doublet L a and the right-handed f a R^aR- Higgs doublets as well as new particles S,a,a',y 
have zero charges. In the limit of exact symmetry, the Higgs doublet and the singlet fermion 

5 can couple only with the electron neutrino, reproducing the matrix (|j). The couplings for 
the second and third generations allowed by U(l) h : 

W = ^TT L ^ + ^h L ^H2 + ^-Ma (19) 
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lead to a mass matrix in {p^v T1 v^pf) basis 



/ m M 

m T 

m M M MT 

V m T M MT J 



(20) 



The mass matrix of charged leptons is diagonal. The diagonalization of (|20|) results in ZKM- 
type (Zeldovich-Mahmoud-Konopinsky) light neutrino formed by ~ and ~ v T components 
with mass 

m 2 = -m 3 = — . (21) 

For ~ 1 GeV, m T ~ 100 GeV and M MT ~ 3 ■ 10 10 GeV one gets m 2 ~ 3 eV which is 
required for the HDM components. In the limit of exact horizontal symmetry v e —S and v^—v^ 
form two unmixed blocks and in particular, m MS = m TS = 0. 

Family symmetry can be conserved at high scale but can be explicitly broken by interactions 
with Higgs doublets. Such breaking could be induced spontaneously also by introducing new 
Higgs doublets with non-zero U(l) h charges (±1 or 2) or by non-renormalizable interactions 
of the type: L e h / ^H 2 cr^/M, where has the charge +1 and acquire the VEV at large scale, 
(o- M ) ~ 10" 4 M. 

Violation of U(l) h leads to mass splitting in v^-v T system as well as to mixing between 
v e —S and v }1 —v T blocks. Consider the phenomenological consequences of introducing U(l) h 
violation separately in different sectors of the model. 

(1) The non-diagonal Dirac mass terms vn^v^v^ + m TM z/ r z/^ + h.c. result in mass-squared 
difference 

Am 2 , * 4m ^ mi . (22) 
For the atmospheric neutrinos one needs Am^ ~ 10~ 2 eV 2 , then for m 2 ~ 2eV and m M ~ 1 



GeV, it follows from (£2|) that m T/Lt should be very small: ~ (0.5 — 1) MeV. Mixing of and 
u T is practically maximal. 

(2) The introduction of a diagonal element in the Majorana sector; e.g., M T v c T v c T , gives 

A < - 2 S (£) < • < 23 > 

and to have Am 2 3 ~ 10~ 2 eV 2 with m il jm T ~ 2 ■ 10~ 2 , one needs M T /M^ T ~ 0.1. 



9 



(3) To get Ve" v n mixing one can introduce the Dirac mass terms m eil p e u'^-\-m t j, e P tl ul + h.c 
Present sensitivity region of KARMEN and LSND: sin 2 29 efl ~ (3 — 5) ■ 10~ 3 corresponds to 
m eA1 /m M ~ 3 • 1CT 2 , and consequently to m eM ~ 30 MeV. In this case v^-S mixing will also be 
generated with tan# MS ~ (m^m^jim^mj) ~ 3 • 10~ 5 which is far below the NS bound. 

(4) Violation of L^l) -symmetry implies in general a non-diagonal mass matrix for the 
charged leptons. In this case the lepton mixing matrix is the product, V — V v ■ V l , where V v 
and Vi diagonalize the mass matrices of neutrinos and charged leptons, respectively. Let us 
suppose for simplicity that the effects of U(l) h violation come from Vi only (V y has two-block 
structure as before), and moreover VJ mixes essentially the first and the second generation 
with the angle Q[. Then the oscillations v e are expected with Am 2 ~ m\ and the 
depth sin 2 26*;. Also mixing between S and appears, so that the z/ M — > S oscillations with 
Am 2 ~ m 2 will have the depth 

sin 2 29 Sfl ~ sin 2 29 es ■ sin 2 9 t . 

For 9 es and 9i fixed by solar neutrino data and the LSND/KARMEN sensitivity, one finds 
sin 2 29 Sfl ~ (1 — 5) ■ 10~ 6 which can satisfy the NS bound. This model realizes the scenario 
described in [|J. 

The v v rv T mass splitting can be generated without explicit U(l) h violation. The modified 
U(l) h charge prescription in the model ( F2" ) with a': (—1,2,0,-2,2,1,0) for {y%,v c , v°,a,cr',S,y) 
allows for the superpotential, 

rxr m e err "V ctj m r CTT 

M e c c I , M P-r c c M T c m es 

+ W) v ' v ' a ~W)"^ T " ~" T " T JH) Ve ' (24) 



For u e -S it reproduces the matrix (^), whereas for v^-Vt system one gets the matrix ( p0|) with 
non-zero M T v^.v° term, thus generating mass splitting fl23"l) . However, the blocks v e -S and 
v^-Vt remain decoupled, and thus no observable effect in KARMEN/LSND is expected. 

5 Conclusion 

We suggest that a light singlet fermion S whose existence is hinted by some neutrino observa- 
tions may have its origin beyond neutrino physics. Such a fermion can however be incorporated 
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into the standard see-saw picture, where interactions of S with the heavy right-handed neutri- 
nos can generate its mixing with the light neutrinos. Such a mixing allows an understanding of 
the lightness of S without ad hoc introduction of very light scale. The mixing mass parameter 
m es ~ (0.02 — 0.3) GeV leads to the mass of the singlet and its mixing with electron neutrino 
in the region mi ~ (2 — 3) • 10~ 3 eV and sin 2 29 es ~ (1 — 15) • 10~ 3 , where the u e — > S resonance 
conversion gives a good fit of all solar neutrino data. 

Super symmetry can provide a framework within which the existence and the desired prop- 
erties of such a light fermion follow naturally. There is a number of models with singlet 
scalars which acquire VEV and are introduced to break symmetries such as lepton number 
and Peccei-Quinn symmetry, or to generate /x-term, etc.. However, identifying S with the 
fermionic superpartner of such scalars implies violation of R-parity, and further complica- 
tion of model. We have considered a specific example with S identified as the majorino. It 
may be possible to suppress the mass of S generated after SUSY breakdown by introducing 
non-minimal Kahler potentials. 

The conservation of R-parity requires for the fermion S to be a component of singlet su- 



perfield which has no VEV. This allows to construct simple model ( fL2|) in which the properties 
(mass and mixing) of S follow from the conservation of R-symmetry. The singlet field is mixed 
with RH neutrinos by the interaction with the field y which can acquire VEV radiatively after 
soft SUSY breaking. The model can naturally incorporate the spontaneous violation of Peccei- 
Quinn symmetry or/and lepton number. The fields involved can spontaneously generate the 
/x-term. 

Approximate horizontal (family) U(l) h symmetry as in ([19]) provides simultaneous ex- 
planations for the predominant coupling of S to the first generation (thus satisfying the NS 
bound) and for the pseudo-Dirac structure of v iir v T needed in solving the atmospheric neu- 
trino and hot dark matter problem. Breaking of U(l) h can be arranged in such a way that 
the parameters of — > v e oscillations are in the region of sensitivity of LSND and KARMEN 
experiments. 

Future solar neutrino experiments will allow to prove or reject the hypothesis of the v e — ► S 
conversion in the Sun |23[ and thus to test the models elaborated in this paper. 



Note added: When our work was practically accomplished we encountered the paper p4 | 



discussing non-supersymmetric model based on discrete symmetry in which sterile neutrino 
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mixes with usual light neutrinos via RH components. Our results have been reported at XXX 
Rencontres de Moriond, March 11-18 (1995), Les-Arcs Savoie, France (to be published). 
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